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Abstract 

One of the important tools to evaluate the ecological quality of surface water is the Macrophytes indi- 
ces based on the bioindication capacity of aquatic plants. In Mediterranean rivers (France, Spain, and 
Portugal), the development of some macrophytes indices like Indice Biologique Macrophytes Riviéres 
(IBMR), the biological metric score (BMS), as well as the Fluvial Macrophyte Index (IMF) are founded 
on the determination of the indicator values of the floristic reference lists. 

The aim of this study was to test the impact of the eco-Mediterranean differences (from one country 
to another) on the indicator taxa by comparing the indicator values of the Euro- Mediterranean mac- 
rophyte indices. With this in mind, we explore the possibility of the introduction of the Euro-Mediter- 
ranean macrophytes-based indices in Morocco (i.e. the hydrological basin of Sebou (HBS)) as a part of a 
preliminary attempt to develop the first Afro-Mediterranean macrophyte index. 

We confirm that the ecological amplitude and species optima vary between Mediterranean ecore- 
gions, and indicator taxa differ between countries: There are medium to small correlations between Medi- 
terranean indices: IBMR/BMS (p = 0.000, R? = 0.57), IMF/BMS (p = 0.000, R? = 0.34), and IBMR/IMF 
(p = 0.000, R’ = 0.30). Five species exhibit major differences in indicator values: Zannichellia palustris 
and Potamogeton pectinatus have more eutrophic indicator values in France (IBMR) than in Spain (IMF). 
Potamogeton nodosus, Amblystegium riparium and Lycopus europaeus have broader ecological amplitudes in 
Portugal (BMS) than in France (IBMR) and in Spain (IMF), where it is restricted to eutrophic conditions. 


Furthermore, the three indicator systems include different indicator-taxon numbers. 
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The comparison of the HBS elaborated list with the Euro-Mediterranean indices revealed the low 
level of common taxa approximately 6.76% of all indicator species used in the French index (IBMR), 
10.48% in the Portuguese index (IMF) and 12.38% in the Spanish index (BMS). 

These results show the inadequacy of the trophic indices approach with the HBS conditions and thus 
the need for the development of an index based on biotic indices approach. 
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Introduction 


Due to their high sensibility to different environmental stresses and their ability to as- 
sess the dynamic and the cumulative effects of different stressing factors, macrophytes 
species are considered good bioindicators. This bioindication power of macrophytes has 
generated a proliferation in the number of macrophyte-based indices in the last decades. 

At the present time, the approaches for estimating macrophyte communities’ qual- 
ity in the Mediterranean rivers are: 


e ‘The approach based on the assumption that environments that have not been im- 
pacted have a greater diversity of species than degraded environments (commu- 
nity structure approaches): Indice di Biodiversita’ Riparia (IBR) (Maggioni et al. 
2009) in Italy is based on biodiversity of macrophytes on the banks. (Patrick 1977) 
proved that assemblages with similar diversity scores could represent streams with 
significantly different chemical conditions. 
e ‘The Biotic indices approach based on the assumption that biological assemblages 
in impaired sites should be different from those in reference sites: 
= ‘The Iberian multimetric plant index (IMPI) (Ferreira et al. 2005), in the Ibe- 
rian Peninsula (Portugal, Spain). 

= ‘The Riparian Vegetation index (RVI) (Aguiar et al. 2009) in Portugal. 

= River Macrophyte Index (RMI) (Kuhar et al. 2011) in Slovenia, based on the 
relative abundance of sensitive and/or tolerant taxa. 

e ‘The approach based on indicator values calculated for an elevated number of 
aquatic species, according to the species’ relative sensitivity and tolerance to nu- 
trients and/or to other abiotic stress factors. The Indices designed to respond to 
nutrient enrichment using indicator species in Mediterranean rivers are: 

© ‘The Indice Biologique Macrophytique en Riviéres (IBMR): developed in 
France by (Haury et al. 2006) for assessing water trophy and organic pol- 
lution and calculated using the following formula: 


n 
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where C’si is the specific rate of trophic level— ranged from 0 (heavy or- 
ganic pollution and heterotrophic taxa) to 20 (oligotrophy); £7 represents 
the coefficient of ecological amplitude: Coefficient 1, representing wide 
amplitude, covered three classes of trophy, and coefficient 3, representing 
a very limited amplitude, was restricted to just one class; Ki is the scale of 
cover, going from +1 t0-5 (1: <0,1%;,2: 0,1 — <1%; 3: 1= = 10%: 4: 10-= 
<50%; 5: >50%). 

© ‘The biological metric scores (BMS): developed by (Dodkins et al. 2012) 
in Portugal. This index is the mean of the species scores that occur at that 
site, weighted by their cover, i.e. the Weighted Averaging (WA) equation 
(Braak and Looman 1986): 


Cixi 
: n 
> Ci 


where S = site score, n = number of species; Ci = cover scale value of spe- 


s= 


cies i; and Qi = score of species i. The cover scale values used to weight the 
mean were: 0 (for 0% macrophyte cover relative to the channel area), 1 
(<1% cover), 2 (<5% cover), 5 (<33% cover) and 6 (>33 cover). 

o ‘The index of macrophytes (IM), the Macroscopic Aquatic Vegetation 
Index (IVAM) and The Fluvial Macrophyte Index (IMF) (Alcaraz et al. 
2006; Flor-Arnau et al. 2015; Suarez et al. 2005) in Spain. The Fluvial 
Macrophyte Index (IMF) is calculated using the following formula: 


)_ Eix Kix Csi 
1M .F =+— 
) Bix Ki 
where £7 is the coating of the taxa at the station -range: 1-5; 1 (<0.1%), 
2 (0.1-1%), 3 (1-10%), 4 (10-50%), 5 (> 50%); Csi is the sensitivity 
value for eutrophy (range: 1—20); £7 is the value of stenoicity or ecological 
amplitude (range: 1-3). The IMF score is obtained from the formula of 
Zelinka and Marvan (1961). 


Taking into consideration that the development of macrophytes assemblages 
strongly depends on a variety of abiotic and biotic factors and it is assumed that the 
most important of them are nutrient concentrations (Dodkins et al. 2012; Robach et 
al. 1996; Schneider et al.2000; Szoszkiewicz et al. 2006; Thiebaut et al. 2002; Whitton 
1975), and hydromorphological characteristics, such as altitude, flow velocity, water 
depth, width of river bed and type of substrate (Balazi and Hrivnak 2017), the overall 
purpose of this paper is to investigate the influence of localized hydromorphological 
differentiation for the bioindication of macrophytes in Mediterranean countries. In 
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particular, we focus on the following question: Is there evidence of a role of hydro- 
morphological differentiation in the diversity of macrophyte taxa included in Mediter- 
ranean indices? Is there any evidence for the impact of ecoregion differentiation on 
the macrophytes indicator values? In other words, are the macrophytes more impacted 
by trophic status or by the hydromorphological characteristics of each Mediterranean 
country? Is there any possibility to adopt and/or adapt any Euro-Mediterranean mac- 


rophytes-based indices in Morocco (HBS)? 


Methods 


All currently used and published Mediterranean macrophyte indices based on species 
indicator values for assessment of river trophic status are included in this study. We 
didn't take into consideration indices with low taxonomic rank resolution (family and 
order): Macroscopic Aquatic Vegetation Index (IVAM) and the index of macrophytes 
(IM). Three macrophyte indices meet the above-indicated criteria: The Fluvial Mac- 
rophyte Index (IMF), the Biological metric scores (S), and Indice Biologique Macro- 
phytes Riviéres (IBMR). 

Comparison of species indicator values between different Mediterranean indices 
was performed using correlation analysis. 

An extensive field survey of macrophytes communities (aquatic and riparian spe- 
cies) in HBS and its tributaries (39 stations) has been carried out. Identification of 
the macrophytes was taken using field identification guides (Ahayoun et al. 2007; 
Coudreuse et al. 2005; Fennane et al. 1999; Fennane et al. 2007; Valdés 2002). 

In order to ensure comparability of species, taxa names were screened for syno- 
nyms and harmonized if necessary. 


Results 


Mediterranean indices comparison 


The most striking results to emerge from Mediterranean indices comparison are: 


IBMR compared to IMF 


A total of 68 species are included in both IBMR and IME Half of these species have an 
IMF value between 16-18 (Figure 3). The indicator values are significantly correlated (p = 
0.000, R*= 0.30) (Figure 2). Two species differ from the regression curve. In the two cases 
the IBMR value is lower than the IMF (Zannichellia palustris, Potamogeton pectinatus). 

A total of 158 taxa have only an IBMR, but not an IMF indicator value, while 56 
taxa have only IMF indicator value but not an IBMR. 
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Figure |. A families B classes, and C orders of macrophytes species included in the Mediterranean 
trophic indices: IBMR, IMF, and BMS. 
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Figure 2. Polynomial regression of A IBMR and IMF B IBMR and BMS C IMF and BMS. 


IBMR compared to BMS 


A total of 47 species are included in both IBMR and BMS. The indicator values are 
significantly correlated (p = 0.000, R’ = 0.57). Two species differ from the regression 
curve. In the two cases the IBMR value is lower than the IMF (Amblystegium riparium, 
Potamogeton nodosus). 

A total of 179 taxa have only an IBMR, but not a BMS indicator value, while 58 
taxa have only an IMF indicator value but not an IBMR. 
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Figure 3. boxplots: indicator values of species that are included in A (IBMR) and (IMF) B (IMF) and 
(BMS) C (IMF) and (BMS). 


IMF compared to BMS 


A total of 35 species are included in both IMF and BMS. The indicator values are 
significantly correlated (p= 0.000, R7=0.34). One species differs from the regression 
curve. In this case the IMF value is lower than the BMS (Lycopus europaeus). 

A total of 89 taxa have only an IMF, but not BMS indicator values, while 70 taxa 
have only BMS indicator values but not IBMR. 


HBS macrophytes compared to European trophic indices 


Our field work and analysis revealed that a limited number (23 indicator species) of 
macrophytes recorded in HBS are utilized as bioindicators in biological monitoring 
for the ecological status assessment in rivers in Euro-Mediterranean countries (Table 
1). Fourteen species are used in IBMR, thirteen species in BMS and IME This limited 
number of indicator species represents only 6.76% of all indicator species used in the 
French index (IBMR), 10.48% in the Portuguese index (IMF) and 12.38% in the 
Spanish index (BMS). 


If we extend our analysis to other European indices ie.: 


e The British index: The Mean Trophic Rank (MTR), there are only ten spe- 
cies of HBS that have MTR indicator value: Berula erecta, Elodea canadensis, 
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Table |. The list of aquatic taxa of HBS that are included in IBMR, BMS, IME 


Species IBMR BMS IMF 

Csi Ei Qi Csi Ei 
Agrostis stolonifera 10 1 12 2 
Arundo donax 1 
Berula erecta 14 wy 
Elodea canadensis 10 2 1 
Epilobium hirsutum 2 4 1 
Equisetum ramosissimum 18 3 
Helosciadium nodiflorum 10 1 a) 4 1 
Hygrohypnum luridum 19 3 
Lemna gibba 5 3 2 8 2 
Ludwigia palustris 5 
Mentha aquatica 12 1 S, 12 2 
Mentha longifolia 18 3 
Mentha pulegium 4 
Nasturtium officinale 11 1 2 8 2 
Phragmites australis 9 2 1 
Potamogeton nodosus 4 3 3 
Potamogeton pectinatus 2 2 8 3 
Ranunculus bulbosus 4 
Rumex conglomeratus 8 2 
Scrophularia auriculata 4 1 
Typha angustifolia 6 2 
Veronica beccabunga 10 1 3 12 3 
Zannichellia palustris 5 1 16 3 


Helosciadium nodiflorum, Hygrohypnum luridum, Lemna gibba, Nasturtium 
officinale, Phragmites australis, Potamogeton pectinatus, Typha angustifolia, 
Zannichellia palustris. 

e The German index: Trophic Index of Macrophytes (TIM), there are only eight 
species of HBS that have TIM indicator value: Berula erecta, Elodea canadensis, 
Mentha aquatica, Nasturtium officinale, Potamogeton nodosus, Potamogeton pectina- 
tus, Veronica beccabunga, Zannichellia palustris. 


All these species are included in the Euro Mediterranean indices, especially in the 
French index. 

One of the most common species used in European countries’ indices (MTR, 
TIM, IBMR, IMF and BMS) and taking place in HBS is Nasturtium officinale. 

Based on IBMR index we have in HBS some species representing wide amplitude 
(Ei = 1): Mentha aquatica, Nasturtium officinale, Agrostis stolonifera, Helosciadium nodi- 
florum, Veronica beccabunga, Zannichellia palustris. And some species representing a 
very limited amplitude (£7 = 3): Hygrohypnum luridum, Lemna gibba, Potamogeton no- 
dosus. Furthermore, some species indicating hypertrophic conditions (e.g. Potamogeton 
pectinatus, Potamogeton nodosus, Csi = 2-4) and others indicating oligotrophic condi- 
tions (e.g. Hygrohypnum luridum, Csi = 19). 

Based on BMS index, species associated with high conductivity and nutrient en- 
richment (Qi = 1) are: Elodea canadensis, Phragmites australis, Arundo donax. 
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IMF index reveals some species representing wide amplitude (Zi = 1): Epilobium hir- 
sutum, Helosciadium nodiflorum, Scrophularia auriculata. Species representing a very lim- 
ited amplitude (Ei = 3): Equisetum ramosissimum, Mentha longifolia, Veronica beccabunga, 
Potamogeton pectinatus, Zannichellia palustris. Some species indicating hypertrophic con- 
ditions (e.g. Epilobium hirsutum, Helosciadium nodiflorum, Scrophularia auriculata, Csi = 
4) and others indicating oligotrophic conditions (e.g. Equisetum ramosissimum, Men- 


tha longifolia, Csi = 18). 


Discussion 


The most obvious difference between the three indicator systems is the number of in- 
cluded indicator taxa: IBMR (226), IMF (124), BMS (105), and TIM (49). 

The IMF and the BMS have the fewest species in common (35 common taxa com- 
pared to 47 between IBMR and BMS and 68 between IBMR and IMF). 

The allocation of the trophic values was based on empirical studies (correlation 
between species occurrence and impact parameters), literature data and expert opinion 
in TIM and IBMR. In BMS and IME, the trophic values were determined only by 
empirical studies. 

IBMR and BMS are moderately correlated (R2=0.57). The worst correlation oc- 
curs between IBMR and IMF (R2=0.30). 

In France (IBMR), Zannichellia palustris and Potamogeton pectinatus have more 
eutrophic indicator values than in Spain (IMF) (Figure 2). Zannichellia palustris is 
commonly associated with nutrient-rich conditions (Vukov et al. 2018) as well as 
Potamogeton pectinatus. For instance, in Germany (Trophic Index of Macrophytes 
(TIM) (Schneider and Melzer 2003)) and Poland (Macrophyte Index for Rivers 
(MIR)) Zannichellia palustris and Potamogeton pectinatus are used as indicator of eu- 
trophic conditions. However, in the UK (Mean Trophic Rank (MTR) (Dawson et al. 
1999)), those species are seen to be tolerant of eutrophication, or cosmopolitan in 
their requirements (Table 2). 

Potamogeton nodosus, Amblystegium riparium and Lycopus europaeus have more oli- 
gotrophic indicator values in Portugal (BMS) than in France (IBMR) and in Spain 
(IMF) (Figure 2). 

In Poland (MIR), Potamogeton nodosus tends to be used to refer to eutrophic condi- 
tions. In Germany (TIM), it is used as an indicator of eutrophic to polytrophic condi- 
tions, which is consistent with the eutrophic BMS, IBMR and IMF indicator values. 
It is therefore likely that Potamogeton nodosus has a broader ecological amplitude. For 
instance, in Zambia (Ihe Zambian Macrophyte Trophic Ranking scheme (ZMTR) 
(Kennedy et al. 2016)), this species is considered as ubiquitous species, occurring 
across from oligotrophic to eutrophic conditions (Table 3). 

Amblystegium riparium is described as tolerant of eutrophication or cosmopolitan in its 
requirements. So, it is therefore likely that this species has a broader ecological amplitude. 
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Table 2. Zannichellia palustris and Potamogeton pectinatus indicator values in MTR, TIM, and MIR. 


Zannichellia palustris Potamogeton pectinatus 
Mean Trophic Rank (MTR) UK STR = 2 tolerant of eutrophication or STR = 1 tolerant of eutrophication or are 
are cosmopolitan in their requirements. cosmopolitan in their requirements. 
Trophic Index of Macrophytes (TIM) IV = 2.93 meso-eutrophic (m-eu) — IV = 2.88 meso-eutrophic (m-eu) — eutrophic 
Germany eutrophic (eu) (eu) 
Macrophyte Index for Rivers (MIR) L = 2 eutrophic L=1 eutrophic 


Poland 


Table 3. Potamogeton nodosus and Amblystegium riparium indicator values in TIM, ZMTR, MTR, 


and MIR. 
Potamogeton nodosus Amblystegium riparium 
Trophic Index of Macrophytes IV=3.1 eutrophic (eu) — eu- 
(TIM) Germany polytrophic (eu-p) 


The Zambian Macrophyte Trophic ZTRSsp=(3 U) ubiquitous species, 
Ranking scheme (ZMTR) Zambia occurring across trophic categories 
from oligotrophic to eutrophic 


Mean Trophic Rank (MTR) UK STR = 1 tolerant of eutrophication or 


are cosmopolitan in their requirements 


Macrophyte Index for Rivers L = 3 eutrophic 
(MIR) Poland 


The apparent weak and moderate correlation and the difference of the included 
taxa and their indicator values from one index to another can be attributed to the hy- 
dromorphological characteristics of the Mediterranean rivers. 

331 species are included in the Euro Mediterranean indices (IBMR, BMS and 
IMF) belonging to 98 families, 66 orders and 24 classes. The most diversified fami- 
lies are: Potamogetonaceae, Cyperaceae, Ranunculaceae, Amblystegiaceae, Typhaceae, 
Plantaginaceae, Characeae, Poaceae, Hydrocharitaceae, Apiaceae, Juncaceae (Figure 
1). The most used genera are: Potamogeton (19 species), Ranunculus (19), Sparganium 
(9), Fissidens (8), Juncus (8), Carex (7), Callitriche (7), Chara (6), Equisetum (5), 
Montia (5) and Najas(5). These indices include some species of Chromista, Bacteria 
and Fungi (Table 4). 

The comparison of the HBS elaborated list with the Euro-Mediterranean indices 
revealed the low level of similarity between HBS community species and the floristic 
reference of the French index (IBMR), the Portuguese index (IMF) and the Spanish 
index (BMS). 

Furthermore, there is a limited number of HBS aquatic species (31 species), which 
is in agreement with previous research (Benamar and Maissour 2014). 

The high level of aquatic species in France and the low-level of aquatic species in 
HBS compared to the Euro-Mediterranean countries can be ascribed to the climate tran- 
sition from thetemperate climate of central Europe to the arid climate of northern Africa 
(Giorgi et al. 2008). These Afro-Mediterranean conditions deeply affect stream flows 
(mixture of perennial and intermittent rivers) and the occurrence of aquatic species. 
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Table 4. List of Chromista, Bacteria and Fungi taxa used in Euro-Mediterranean indices. 


kingdom species IMF IBMR 
Bacteria Nostoc + + 
Oscillatoria + + 
Phormidium + + 
Sphaerotilus tt 
Chromista Cymbella + 
Leptomitus 4 
Melosira + Z 
Tribonema + + 
Vaucheria + + 
Fungi Collema dichotomum + 
Dermatocarpon luridum + 


These results demonstrate the inadequacy of the trophic indices approach espe- 
cially with the HBS conditions and in general in the Afro-Mediterranean region, and 
thus the need for the development of an index based on biotic indices approach taking 
into consideration also the riparian species. 

The Biotic indices approach, which is originally developed by Karr and Dudley 
(1981), is a widely used method for evaluating anthropogenic pressures on aquatic and 
wetland ecosystems: Floristic Quality Assessment Index (FQAI) (Lopez and Siobhan 
Fennessy 2002), Integrity Biotic Index (IBI) (Miller et al. 2006), Iberian Multi metric 
Plant Index (IMPI) (Ferreira et al. 2005), Index of Plant Community Integrity (IPCI) 
(DeKeyser et al. 2003), Index of biotic integrity in Itanhaém (MIBI-ITA) (Umetsu et 
al. 2018), Plant Index of Biotic Integrity (PIBI) (Simon et al. 2001), Plant-based index 
of biotic integrity (PIBI) (PIBI(M)) (Moges et al.2016), Riparian Forest Quality index 
(QBR) (Munné et al. 2003), Riparian Quality Index (RQI) (Del Tanago et al. 2006; 
Gonzalez del Tanago and Garcia de Jalon 2006), Vegetation Index of Biotic Integrity 
(VIBI) (Mack 2007), and Vegetation-based index of biotic integrity (VIBI(Y)) (Yang 
et al. 2018). 

Among the potential characteristics of the aquatic vegetation (candidate met- 
rics) that can be responsive to disturbance in HBS are: diversity, species habitat, 
life cycle, life form, nutritional resources, riparian structure, and species tolerance 
(Table 5). 

Future work will involve the selection of the reference sites. This is because the 
reference sites provide the baseline information to detect the deviation of a metric 
from a natural or least-disturbed condition. And the selection of suitable metrics in our 
context. So, we need to evaluate the ability of every potential candidate metric in terms 
of its ability to distinguish reference (undisturbed or least-disturbed) from impaired 
(moderately or heavily disturbed) sites. Only the metrics showing significant difference 
between reference and impaired sites will be selected as the IBI-HBS metrics (Yang et 
al. 2018). The next step is to score the selected core metrics. 
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Table 5. Potential candidate and core metrics in IBI- HBS. 


Candidate metrics Expected response to —_— River and wetland indices 
decreasing quality 
Diversity 
Species richness Decrease FQAI, PIBI, IPCI, MIBI-ITA 
Species habitat 
% Endemic species Decrease FQAI 
% Native species Increase PIBI(M) 
% Exotic species Increase PIBI(M) 
Life cycle 
% Annual species Decrease IMPI, IBI, VIBI, 
% Perennial species Increase VIBI, PIBI, IPCI, VIBI(Y) 
Life form 
% Terrestrial species 
% Hygrophyte species RQI, PIBI 
% Helophyte (emergent species)+ hydrophyte species Decrease IMPI, VIBI, PIBI; MIBI-ITA 


(floating-leaved, free-floating, and submerged species) 
Nutritional resources 


% Ruderal species Increase IMPI 
% Nitrophyllous species Increase IMPI, RQI 
Riparian structure 
% Woody species richness (trees, shrubs, woody climbers) Variable IMPI, IBI, RQI, PIBI(M) 
Species tolerance 
Tolerant species richness Increase PIBI(M), VIBI(Y) 
Sensitive species richness Decrease PIBI(M) 
Conclusion 


We have confirmed that the ecological amplitude and species optima vary between 
Mediterranean ecoregions, and that indicator taxa differ between countries. 

It was found that the trophic indices of the Euro Mediterranean rivers can't be applied 
easily to the Afro- Mediterranean rivers, particularly in Morocco (HBS), and we don’t have 
a good opportunity to enrich the list of indicative species due to the limited number of 
species recognized as bioindicators (23 species) and the limited number of aquatic species. 
So, it seems more appropriate to develop an index based on a biotic-integrity approach. 
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